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ABSTRACT 

Hundreds of substellar companions to solar-type stars will be discovered with the Kepler satellite. Kepler's, 
extreme photometric precision gives access to low-amplitude stellar variability contributed by a variety of 
physical processes. We discuss in detail the periodic flux modulations arising from the tidal force on the star 
due to a substellar companion. An analytic expression for the variability is derived in the equilibrium-tide 
approximation. We demonstrate analytically and through numerical solutions of the linear, nonadiabatic stellar 
oscillation equations that the equilibrium-tide formula works extremely well for stars of mass <1.4M Q with 
thick surface convection zones. More massive stars with largely radiative envelopes do not conform to the 
equilibrium-tide approximation and can exhibit flux variations >10 times larger than naive estimates. Over the 
full range of stellar masses considered, we treat the oscillatory response of the convection zone by adapting a 
prescription that A. J. Brickhill developed for pulsating white dwarfs. Compared to other sources of periodic 
variability, the ellipsoidal lightcurve has a distinct dependence on time and system parameters. We suggest that 
ellipsoidal oscillations induced by giant planets may be detectable from as many as ~100 of the 10 Kepler 
target stars. For the subset of these stars that show transits and have radial-velocity measurements, all system 
parameters are well constrained, and measurement of ellipsoidal variation provides a consistency check, as 
well as a test of the theory of forced stellar oscillations in a challenging regime. 
Subject headings: planetary systems — stars: oscillations — techniques: photometric 



1. INTRODUCTION 

The upcoming Kepler' satellite will continuously monitor 
~10 5 main-sequence stars of mass ~0.5-l.5M© over 4-6 
years with fractional photometric precisions of ~10~ 5 . Such 
high sensitivity, which is unattainable from the ground, will 
allow for the robust detection of Earth-size planets that transit 
their host stars, and the measurement of a steroseismic oscilla- 
tions as a probe of stellar structure (e.g.. lBorucki et al"1l2004t 
iBasri et al.l2005l) . These missions will also discover hundreds 
of "hot Jupiters" with orbital periods of < 10 days, revealed 
by th eir transits or reflected starlight (e.g., iJenkins & Dovld 
2003). Continuous observations of these systems are likely to 
show a myriad of novel physica l effects, including D oppler 
flux variability of the host stars (lLoeb & Gaudil l2003). pho- 
tometric dips due to moons or rings around t he planets 
dSartoretfi & Schneider! 11999b iBrown et al l 1200 lb . and the 
impact of additional perturbing p lanets on transit timing 
( Miral da-Escudl l2002t IXgol et al J 120051: iHolman & Murray! 
120051) . The same ideas apply if the companion is a more mas- 
sive brown dwarf, but these are rarely found in close or bits 
around solar-type stars (e.g., Grether & Lineweaver 2006). 

Here we scrutinize another mechanism for generating peri- 
odic variability of a star closely orbited by a giant planet or 
brown dwarf. A star subject to the tidal gravity of a binary 
companion has a nonspherical shape and surface-brightness 
distribution. In the simplest approximation, the stellar surface 
is a prolate ellipsoid with its long axis on the line connecting 
the two objects. As the tidal bulge tracks the orbital motion, 
differing amounts of light reach the observer. For a solar-type 
star orbited by a perturbing companion of mass M p with pe- 
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riod P or b, the expected fractional amplitude of this ellipsoidal 
variability is ~10" 2 (M / ,/M Q )(l day/f orb ) 2 . This effect has a 
long his tory in the stu dy of eclipsing binary stars (see the re- 
view by lWilsorJ l994). but was mentioned only recently in the 
e xoplanet context. 

lUdalski et al.ld2002l) . lDrakel(l2003l) . and lSirko & Paczvnskil 
(2003) noted that if ellipsoidal light variations are detected 
from the ground, where the fractional photometric precision 
is >10~ 3 , then the perturber must be fairly massive (e.g., 
>0.1M Q ). They offered this idea as a test to distinguish be- 
tween planetary transits and eclipses by low-mass stars. The 
superior sensitivity of Kepler offers the possibility of mea- 
suring ellipsoidal variability induced by giant planets (M p ~ 
1 0~ 3 -10~ 2 My) w ith orbital periods of <10 days. 

iLoeb&Gaudl (120031) compare the ellipsoidal variability in- 
duced by a planetary companion to flux modulations arising 
from reflected starlight and the Doppler effect. The three am- 
plitudes are similar when the companion has an orbital period 
of <3 days and an optical albedo of <0.1. In a sufficiently 
long observation it should be possible to separately extract 
each of the signals, since their Fourier decompositions are dis- 
tinct. Precise physical modeling of the ellipsoidal lightcurve 
could provide an independent constraint on the mass of the 
companion, as well as important clues regarding stellar tidal 
interactions. 

Ellipsoidal variability is typically modeled under the as- 
sumption that the distorted star maintains hydrostatic balance 
and precisely fills a level surface of an appropriate poten- 
tial (e.g., the Roche potential). The measured flux is then 
just an integral of the intensity over the visible stellar sur- 
face, where the intensity includ es the effects of limb darken- 
ing and gravity darkening (e.g., Kopal 1942). This approach 
is strictly valid only when the orbit is circular and the star 
rotates at the orbital frequency, so that a stationary configura- 
tion exists in the coorbital frame. These conditions may not 
be satisfied when the companion has a low mass or long pe- 
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riod, because of the weak tidal interaction. In fact, a state of 
tidal equilibrium may not be attainable i n the case of a plan- 
etary companion (e.g. jRasio et alJ ll996). Equilibrium mod- 
els of ellipsoidal lightcurves do have a realm of validity for 
noncircular orbits and asynchronously rototating stars, and 
have bee n applied successfully to somewhat eccentric bina- 
ries (e.g., ISoszvnskl et al. 2004). However, by construction, 
such models ignore fluid inertia and the possibility exciting 
normal modes of oscillation, effects that may be of critical 
importance in a wide range of observationally relevant cir- 
cumstances. Here we apply the machinery of linear stellar 
oscillation theory to the weak tidal forcing of stars by sub- 
stellar companions. Conceptually, our investigation bridges 
Kepler's planetary and astroseismology programs. 

Section 2 describes the geometry of the problem, provides 
quantitative measures for the strength of the tidal interaction, 
discusses our simplifying assumptions, and presents the math- 
ematical framework for calculating ellipsoidal variability. In 
§ [3] we consider the equilibrium-tide approximation and de- 
rive an analytic expression for the ellipsoidal lightcurve. A 
brief review of von Zeipel's theorem and its limitations is 
given in § |U Tidally forced, nonadiabatic stellar oscillations 
are addressed in § [5] where we argue for a simple treatment 
of perturbed surface convection zones, use this prescription to 
calculate the ellipsoidal variability of deeply convective stars, 
estimate analytically the surface flux perturbation in mainly 
radiative stars, and show select numerical results. Our main 
conclusions are summarized in § [6j We conclude in §[7] with 
remarks on the measurement of ellipsoidal oscillations in the 
presence of other sources of periodic variability. 



2. PRELIMINARIES 

Consider a star of mass M and radius R is orbited by a sub- 
stellar companion of mass M p and radius R p . We work in 
spherical coordinates (r, 8, </>) with the origin at the star's cen- 
ter and the pole direction (9 = 0) parallel to the orbital angular 
momentum vector. The orbit is then described by (d, ir/2, (j> p ), 
where d and <f> p are, respectively, the time-dependent or- 
bital separation and true anomaly; <fi p = marks the phase 
of periastron. We assume that the orbit is strictly Keple- 
rian with fixed semimajor axis a and eccentricity e, such that 
d = a(l-e 2 )/(l +eco&(f> p ). The direction to the observer from 
the center of the star is (9„, </>„), so that the conventional orbital 
inclination is / = %— 6 a . 

We imagine that the gravity of the companion raises nearly 
symmetrical tidal bulges on opposite sides of the star that ro- 
tate at the orbital frequency. A rough measure of both the 
height of the tides relative to the unperturbed stellar radius 
and the fractional amplitude of the ellipsoidal variability is 
given by the ratio of the tidal acceleration to the star's surface 
gravity: 
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where Mj ~ 10 3 M Q is the mass of Jupiter, P„ = 
2w(R 3 jGMfl 1 = 2.8[(/?//? ) 3 (M Q /M)] I / 2 hr is the dynami- 
cal time of the star. For main-sequence stars with R/Rq ~ 
M/Mq, we see that e oc M p MP~^. The maximum value of 
£ is attained when the companion fills its Roche lobe at an 
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where we have applied a fixed value of R p = 0.1 Rq, appro- 
priate for both giant planets and old brown dwarfs. Note that 
£max ~ 1 for massive brown dwarfs (M p /Mj ~ 80). Hereafter, 
we consider only cases with £<1. 

For orbital periods as short as ~ 1 day, tidal torques on the 
star from a planetary companio n are rather ineffec tive at alter- 
ing the stellar rotation rate (e.g jRasio et al.ll9 96). Therefore, 
as already mentioned in § Q] we should not generally expect 
the star to rotate synchronously with the orbit, and so there 
is no frame in which the star appears static. This holds when 
the orbit is circular, and is obviously true when the there is a 
finite eccentricity. In fact, ~30% of the known exoplanets 4 
with P 01 b < 10 days have eccentricities of >0.1. Small vari- 
able distortions of the star from its equilibrium state, due to 
a combination of asynchronous rotation and orbital eccentric- 
ity, should be viewed as waves excited by the tidal force of 
the companion. Our task is to study such tidally forced stel- 
lar oscillations in the linear domain in order to understand the 
corresponding lightcurves. 

In order to greatly simplify the mathematical description of 
the stellar oscillations, we assume that the star is nonrotating 
in the inertial frame. When the stellar rotation frequency is 
nonzero, but much smaller than the tidal forcing frequency, 
the effect of rotation is to introduce fine structure into the os- 
cillation frequency spectrum, and cause the oscillation eigen- 
functions to be slightly modi fied as a result of the Coriolis 
force (for a discussion, see lUnno et al.|[l989h . Tidal pump- 
ing of a slowly rotating star by an orbiting companion has a 
dominant period of P or b/2 — a few days in the cases of inter- 
est. By contrast, single solar-type stars with ages > 1 Gyr tend 
to have rotation period s of >10 days (e.g., Sku manic 



jyr tend 



iPace & Pas quini 2004); the Sun has an equatorial rotation pe- 
riod of ~25 days. Slowly rotating stars with masses of ~1 M Q 
are prime targets for Kepler, since they exhibit low intrinsic 
variability. Based on this selection effect, and the inability of 
tidal torques to spin up the star, our assumption of vanishing 
stellar rotation seems generally justified. 

The general framework for calculating the measurable flux 
modulations associated with ellipsoidal stellar oscillations is 
as follows. We consider small perturbations to a spherical, 
nonrotating background stellar model, such that fluid ele- 
ments at equilibrium position x are displaced in a Lagrangian 
fashion to position x+£. Variations in the measured flux from 
an oscillati ng star arise from tw o physically distinct contribu- 
tions (e.g., Dziembowski 1977h : (1) changes in the shape of 
the star due to radial fluid displacements = £ • e r , where e r 
is the radial unit vector, and (2) hot and cold spots generated 
by local Lagrangian perturbations AF to the heat flux. Our 
main task in §§[3]and|5]is to compute £ r and AF according to 
the relevant physics. 

Given the dependences of £ r and AF on (r,6,<fi), it is 
straightforward to compute the time varying component of the 
measured flux. The flux 5 received from a star at distance D is 

4 |http : //vo . obspm. f r/exoplanetes/encyclo/encycl . html| 
5 Our calculations concern the bolometric flux, although is relatively 
straightforward to modify the analysis for narrow-band measurements. 
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(e.g., iRobinson etalll 19821) 
1 



& = ^2 / dSn-n a F h(n ■ n e ) 



(3) 



where dS is an area element at the stellar photosphere, F is 
the net flux of radiation out of the surface element, h is the 
limb-darkening function, n and n„ are unit vectors normal to 
the surface and toward the observer, respectively, and the in- 
tegration is over the visible stellar disk. Vertical displacement 
at the surface yields chang es in & through chang es in surface 
area and n • n . Following Dziembowskj (1 19771) . we expand 
£ r and AF in spherical harmonics, 



(4) 
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and carry out the appropriate linear expansions to obtain the 
fractional variability 
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Here and AF£ are components evaluated at the surface 
(r = /?) and in the direction of the observer: 
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The terms £^ and q are given by 



d^JL^iPlh , Q : 



dji{\- ji )- 



h+fj,- 



dh 



(7) 
(8) 

,(9) 



/o jo ^A* V ^A*, 

where /i = n • n„, the are ordinary Legendre polynomials, 

and h(fi) is normalized such that L d/ifih = 1. The linear 
limb-darkening function is 

/i(/i)= -^-[1- 7 (1- M )] ; (10) 
(3-7) 



more general nonlinear functions of fj, (e.g., IClaretl 120001) 
will not be considered here. The classical Eddingto n limb- 
darke ning function is h = 1 + 3/i/2 (7 = 3/5; e.g., Mihalas 
1970). Table Q] shows shows functional forms and particular 
values of be and a for £ = 2 and 3. 

3. EQUILIBRIUM TIDE 

Vertical displacement of the stellar surface is often accu- 
rately modeled by assuming that the tidally perturbed fluid 
remains in hydrostatic balance. The cause and magnitude of 
the surface flux perturbation is a more complicated affair. In 
this section, we apply a simple parameterization of the flux 
perturbation and obtain a complete set of formulae for com- 
puting the ellipsoidal lightcurve. Subsequent sections provide 
more detailed calculations. In particular, we show in § 15.21 
that stars with deep convective envelopes (the majority of Ke- 
pler targets) have surface flux variations that conform to the 
equilibrium-tide approximation. 



TABLE 1 
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When the tidal forces on the stellar fluid change sufficiently 
slowly, the star can stay very nearly in hydrostatic equilib- 
rium. If the net acceleration required to balance the pressure 
gradient is derivable from a potential, then equilibrium im- 
plies that a fluid element remains on an equipotential surface. 
Since we neglect stellar rotation, there is no centrifugal force, 
and the total potential is the sum of the gravitational potential 
99 from the spherical background stellar model and the per- 
turbing tidal potential U ~ eip -C f. For our analytic work, 
we neglect the modification of ip due to the tide. In general, 
the Eulerian variation Sep should be added to U, as we do in 
our numerical models (see § 15 .41 and the Appendix); we find 
that \8tp\/\U\ ~ 10" 2 . 

In the absence of tidal forces, a given fluid element sits at 
equilibrium position x with total potential ip(x). Gentle inclu- 
sion of the tidal potential causes the fluid element to move to 
position x+£ while preserving the value of the total potential. 
This is expressed mathematically by 

<p(x) = <p(x+€) + U(x+$,t) 

= tp(x) + Z-V<p + U(x,t) + 0(e^U) . (11) 

We see that £ • V<p = £ r g, where g = GM r / r 2 is the background 
gravitational acceleration at mass coordinate M r . To first or- 
der, the radial displacement of the equilibrium tide is (see also 
iGoldreich & Nicholson 1989) 



£ r (x,t)~-U(x,t)/g. 



(12) 



which tells us the geometry of the star as a function of time. 
The tidal potential within the star can be expanded as 



GM„ 



(-)Pi(cosVO, (13) 



U(r,6,cj>,t) = -^-£ 

d t - J \d 

e=2 

where cos tjj = sin#cos(</) ;) - 0). There is no £ = 1 term, since 
this would give the acceleration of the star's center of mass, 
which is already incorporated into the orbital dynamics. The 
angular expansion of follows immediately from eq. (fT2l : 

~ r = WXU PeiCOS ^- (14) 

In order to express U and £ r in spherical harmonics, we utilize 
the addition theorem, 



Pe(cosijj) = 



4-TT >. - 

2i+\ 2-", 



Yl m (n/2,<fi p )Yt m (6,<P), (15) 



where "*" denotes the complex conjugate. Note that 
^m(7r/2, 4>p) is nonzero only when t—m is even. For the dom- 
inant 1=2 components of U and £ r , the surface values of U/(p 
and £ r /R are ~e, as expected. 

From eqs. ( fl4l >. and (15[ , the components £°#/R of the 
surface radial displacement toward the observer are immedi- 
ately apparent. As we will see in § [5] the computation of 
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AF/F is, in general, rather technical. However, in the special 
case where the stellar fluid responds adiabatically to a slowly 
varying tidal potential, AFi/F varies in phase with and in pro- 
portion to £ r j jr in the linear approximation of the equilibrium 
tide. Making this assumption, we write AFg/F = -Xg^i/R at 
the surface, where the Xj ar e real constants that depend on the 
stellar structure (see § 15. 11 1. We will see in § |4]that Xi = 1+2 
is a good first guess for radiative stars, and so we might gen- 
erally expect X( to be positive and O(l). 

We now have the ingredients for the fractional variability 
(eq. J6)), and we obtain 



5& ^ (R 

^ 1=2 



ft Pg(COSl/)„) 



(16) 



where fa = (2-Xe)b£-ce, and cosi[j = smd cos(c/) p -(j) ). The 
i = 2 and 3 Legendre polynomials can be expanded as 

P 2 (cosV^)=i[-(3cos 2 /-l) 

+ 3sin 2 / cos 2(cf> p -</>„)] , (17) 

1 r , 
Piicosipo) = - sin/ [-3(5 cos /- l)cos(tf> p - cj) ) 
8 

+ 5sin 2 / cos 3 (0 ; , -</>„)] , (18) 

where we have substituted d„ = it - 1. The Eddington limb- 
darkening formula gives (see Table [U 
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(19) 



It is important to note that fa < when A2 > (see below). 

In eq. ( TToT ). the orbital dynamics are described by the evolu- 
tion of d and <f> p (see §|2|. For a circular orbit, we have d = a 
and (j) p = fit, where fl = 27r/P orb , and t is the time since perias- 
tron (modulo P or b). Example lightcurves with e = 0, 7 = 3/5, 
Xg = 0, and /= tt/2 are shown in Fig.[T]for a/R = {2,4,8, 16}. 
When /J/aCl, the 1 = 2 piece of 8^/^ is a good approxi- 
mation, and the temporal flux variation approaches a pure co- 
sine with angular frequency 20, (see eq. ffTTl ). Because fa < 0, 
the dominant 1 = 2 component of the ellipsoidal variability has 
minimum light when tidal bulge is aligned with the direction 
to the observer. As R/a increases, so does the importance 
of £ > 2 terms and their extra harmonic content, as seen in 
eq. ([H and Fig. [1] 

Additional harmonics in also result from a finite 

eccentricity. At the 0(e) level, signals with frequencies £1 and 
3ft, and amplitudes of ~ee, are present in the 1 = 2 component 
of <5JF/J^, which compete with the I = 3 piece when e ~R/a. 
Notice that when e > the flux is variable even when the orbit 
is viewed face-on (/ = or it), by virtue of changes in d~ 3 = 1 + 
3ecosf2f + 0(e 2 ). For / = 0, we see that /^(cosV'o) vanishes, 
leaving the largest contribution 8& ~ — \.5eefacos(Q,i). 

4. AN ASIDE ON VON ZEIPEL'S THEOREM 

Our equilibrium calculation in the last section used the sim- 
ple prescription AFi/F = -Xi£ r ,i/R- There remains the ques- 
tion of what physics determines AF/F. A common practice 
in empirical studies of close eclipsing binaries — systems that 
tend to be nearly in t i dal eq uilibrium — is to use some vari- 
ant of the Ivon Zeipefl (119241) theorem, which was originally 
formulated for purely radiative, strictly hydrostatic stars. In 
equilibrium, all the thermodynamic variables depend only on 




FIG. 1. — Disk-averaged flux variation (eq. 1161 ) for an edge-on circu- 
lar orbit under the equilibrium-tide approximation (eqs. (6) and 1141 ) with 
AF/F = at the surface. The four curves correspond to a/R = 2 (black), 4 
(red), 8 (blue) and 16 (green). In order to compare the shapes of the curves, 
<5J7Jf has been multiplied by (a/R) 3 (M/M p ). As a/R increases, higher har- 
monics decrease in strength and the lightcurve approaches a pure cosine with 
frequency 2/P or b- The tidal bulge closest to the companion points toward the 
observer at integer values of f/P or b. 

the local value of the t otal potential j>. Thus, the radiative flux 
can be written as (e.g., Hansen & Kawaler 1994) 
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(20) 



where p is the mass density, T is the effective temperature, 
and k(p,T) is the opacity. Equation ( f20b is the essence of 
von Zeipel's theorem, which says that the magnitude F of the 
radiative flux is proportional to the magnitude of the net ac- 
celeration A = I V<&|. When <i> = ip + U (see § [3]), we obtain 
A = g + dU I 'dr+ 0(£ 2 ), so that the Lagrangian flux perturba- 
tion about equilibrium is 



AF 
~F~ 



AA 



Ag 1 dU 
R X dr 



(21) 



where Ag/g = -2£, /r, due to the change in radius at approxi- 
mately constant enclosed mass. Substituting the equilibrium- 
tide result U = -£ r g into eq. d2"Tl ). we obtain the compact ex- 
pression AF/F = -d^ r /dr. Using eq. (fl4l i. we find 



AF e 
~F~ 



1+2) 



(22) 



from which we identify Xe=£ + 2. 

Although the application of von Zeipel's theorem is instruc- 
tive, the underlying physical assumptions are inaccurate for 
slowly rotating main-sequence stars of mass l.O-1.6M with 
tidal forcing periods of days. We are now led to investigate the 
general problem of forced nonadiabatic stellar oscillations. 

5. FORCED NONADIABATIC OSCILLATIONS 

The equilibrium analysis ignores fluid inertia and the exci- 
tation of the star's natural oscillation modes. While this as- 
sumption may be valid near the surface of the star, it does 
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lg Pressure[dyne crrr 2 ] 

FIG. 2. — Important oscillations frequencies and time scales as a func- 
tion of pressure for a 1 Mq main-sequence star. The four curves show the 
Brunt- Vaisalla frequency N (black), Lamb frequency Si (red; £ = 2 is shown), 
inverse thermal time C 1 (blue), and inverse eddy turnover time C? (green). 
Large, real values of N occur in the radiative core and very near the photo- 
sphere, while N < in the convective envelope. Gravity waves propagate 
only where the angular frequency is below both N and Li . The two horizontal 
lines delimit the range of tidal forcing frequencies of interest here. 

not hold deeper in the interior. Gravity waves (g-modes; re- 
stored by buoyancy) can propagate in the radiative interiors 
of Sun-like stars with a range of oscillation periods that in- 
cludes the tidal forcing periods of interest (<3 days). Tidal 
forcing of radiative regions may produce substantial devia- 
tions from hydrostatic balance, as well as large surface am- 
plitudes of AF/F, in particular if resonant oscillations are 
excited. This is especially relevant for main-sequence stars of 
mass M > 1.4— 1.5M© with mainly radiative envelopes. Less 
massive stars (M < 1.3-1.4 M©) have rather deep convective 
envelopes that can block information about the dynamic inte- 
rior from being conveyed to the surface. Here we investigate 
each of these regimes with both analytic estimates and numer- 
ical models of oscillating stars. 

Our calculations employ realistic models of 0.9-1.6M Q 
main-sequ ence stars, co nstructed with the EZ stellar evolu- 
tion code {Paxton 2004), a distilled and rewritten version of 
the program originally created by Peter Eggleton. We adopt 
Solar metallicity and a convective mixing length of 1 .6 times 
the pressure scale height. All stars are evolved to an age when 
the core hydrogen abundance has the Solar value ofXy = 0.35. 
Models with 199 radial grid points are interpolated to yield 
>10 4 points in which the g-mode radial wavelength is well 
resolved in the core. 

Figures |2] and [3] illustrate some of the differences between 
1M© and 1.6M Q stars, and serve to introduce several impor- 
tant physical quantities used in the remainder of this section. 
The Lamb frequency, 



S e = [£(£+!)] 



1/2 £i 



(23) 



is the inverse of the horizontal sound-crossing time scale, 
where c s is the sound speed, and [£(£+ l)]'/ 2 /r = kh is the 
horizontal wavenumber of the oscillation. For fixed chemical 
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FIG. 3. — Same as Fig. [2] but for a 1.6A/q main-sequence star. Note two 
geometrically thin, relatively inefficient (f t h ~ ?ed) convection zones near the 
surface. The spike in N near the center is at the edge of the convective core, 
and signals a steep gradient in the mean molecular weight. 

composition, the squared Brunt-Vaisalla frequency is 

jV 2 ~-J-(V ad -V), (24) 
tip 

where H p = -(dlnp/dr)~ l = p/(pg) is the pressure scale 
height, and V = d In T /d lnp is the temperature gradient 6 (V a( j 
is the adiabatic value). Radiative regions have V a d - V > 
(Af 2 > 0), and Af represents the frequency of buoyancy oscilla- 
tions. In convection zones, V a d-V < andA^ 2 < 0, indicating 
that g-modes are evanescent. When A^ 2 < 0, the time scale 



fed ~ \N\- 



(25) 



approximates the turnover time of convective motions (for de- 
tails and modifications for radiative losses, see, e.g., Kippen- 
hahn & Weigert 1990). A shell of radius r, thickness H p (size 
of the largest convective eddies), and radiative luminosity L 
cools on the thermal time scale 



fth = 



4Trr 2 H p pC p T 



(26) 



where C p is the specific heat at constant pressure. 

The 1 Mq model (Fig.O has one deep convection zone with 
fed <C fth over most of the region, indicating that convection 
very efficiently transports energy and causes the zone to be 
essentially isentropic. By consrast, the 1.6M star (Fig. [3} 
has two thin surface convection zones with f e d ~ fth, and thus 
the radiative and convective fluxes are comparable. Gravity 
waves with frequency u> propagate only in radiative regions 
where u> < N and u> < S^. For the 1M Q star, heat and en- 
tropy generated by g-modes in the radiative interior may be 
strongly mitigated owing to the long thermal time at the base 
of the deep convection zone. On the other hand, g-modes in 
a 1.6M Q star can propagate very near the surface, producing 
qualitatively different results. 

6 Do not confuse the temperature gradient V with the spatial gradient V 
used in §f5] 
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age[Gyr] 

FIG. 4. — Eddy turnover time at base of the convective envelope versus 
stellar age for a range of stellar masses. 

We now go on to elucidate the physics of the flux per- 
turbations. All the analytic and numerical work that fol- 
lows assumes that the tidal potential has the generic form 
U oc rYt m (6 ,<fi)exp(-iojt) with forcing frequency u>. 

5.1. Heat Transfer in a Convective Envelope 

Calculation of the perturbed convective flux in oscillating 
stars is a thorny issue. For the purposes of our study, we argue 
for an especially simple treatment that draws from previous 
work on this subject. Specifically, we modify the prescription 
of Brickhill (1983, 1990; see also Goldreich & Wu 1999a,b), 
which was originally applied to white-dwarf pulsations, into 
a form appropriate for the tidal flow problem. 

In the mixing-length theory of convection, heat is trans- 
ported by eddies with a spectrum of sizes / < H p , speeds vj, 
and turnover times f e d(0 = // v; . The Kolmogorov scalings for 
turbulent motions give v/ oc /'/Vedoc/^and an energy den- 
sity per unit mixing length interval oc/" -1 / 3 . We see that in the 
unperturbed star most of the convective energy flux (ocv 3 at 
scale I) is carried by the largest eddies (I ~ H p ). Convection 
efficiently transports energy when the radiative thermal time 
scale associated with the dominant eddies is much longer than 
t e d. Alternatively, efficient convection implies that the gradi- 
ent of the specific entropy s is small; i.e., dlns/dlnp <C 1. If 
all the convective energy flux F is carried by eddies with mix- 
ing length /, the flux and entro py gradient are related by (e.g., 
Kippenhahn & Weiger3 ll990b 



1 ds 
C p dlnp 



(Z) = (V-V ad ) 



F 

PCs 



2/3 



(27) 



Efficient convection enforces V - V a d *C 1 , which implies 
dlns/dlnp <C 1, since s > C p in the convective regions of 
our background models. 

Gravity waves with the tidal forcing frequency lu are ex- 
cited in the radiative region below the convection zone. Con- 
vective eddies can tr ansport heat during a forcing period only 
if fed < 2it/lu (e.g.. lBrickhilil[T990l: iGoldreich & Wul ll999b). 




age[Gyr] 



FIG. 5 . — Thermal time at the base of the convection envelope versus stellar 
age for a range of stellar masses. 

Inspection of Fig. [2] shows that in the IM Q model, the 
largest eddies have f e d — 20(p//?t>cz) ' 5 days, where pbcz — 
10 13 - 5 dyne cm" 2 is the pressure at the base of convection 
zone. Using the Kolmogorov scaling, the "resonant" length 

/"res for Which CJf e d/27T = 1 is 

3/2 / „ \ 3/4 

1(r MT-E) l^J > (28) 

>1 for 



k 



which is 



2tt/u> 
1 day 
all periods 



2tt/l 



Pba 
P 

' > 



1 day when p < 



10 12 dyne cm , which still encompasses much of the convec- 
tion zone. Now imagine the situation where all the convective 



flux is carried by eddies of size </„ 
for this range of mixing lengths is 



The entropy gradient 



1 ds 



-(/"res) -10" 



/ 2tt/o 



\ 1 day 



P 

Pbc; 



(29) 



jp d In p 

where we have adopted F / pc s ~ 10~ 8 at the base of the con- 
vection zone, as indicated by our 1 M stellar model. 

These arguments suggest that convection is efficient in a 
1 Mq star at the forcing periods of interest even if small "res- 
onant" eddies carry all the energy flux near the base of the 
convection zone. At larger radii, but not too near the photo- 
sphere, convection is both efficient and rapid ((x>f e d/27r < 1) 
over the full spectrum of eddies. Rapid convection on all 
scales I < Hp enforces isentropy in the convection zone, such 
that s an d its Lagrangian pertur bation As are nearly constant, 
as in the iBrickhilH d!983l 11990b picture. While convection at 
the base is rapid only on small scales, it is still highly efficient, 
which yields s ~ constant and further indic ates that As/C p is 
small in magnitude, as we demonstrate in § 15.21 

As the stellar mass increases, the convection zone thins and 
f e d at the base decreases (see Figs.[3]and|4]i. Rapid convection 
holds over the bulk of the convection zone for masses >IMq. 
However, the assumption that the convection is efficient starts 
to break down at 1.4— 1.5 Mq, since f t h ~ f e d at the base (see 
Figs. [3] and |5}. For the full range of stellar masses consid- 
ered here, we assume that s and As are constant in convection 
zones. 
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5.2. Analytic Result for Thick Convection Zones 

In a fully convective star, the emergent luminosity is de- 
termined entirely by the surface boundary conditions. Under 
our assumption that As is constant in the convection zone, the 
perturbed luminosity is likewise a function only of the bound- 
ary conditions. Stars of mass <1.3-1.4M Q have long thermal 
times (w?th 3> 1) at the top of the interior radiative region (see 
Fig. [5J, so that the flux perturbation AF is approximately the 
"quasi-adiabatic" value, derived by ignoring As <x (ujt±)~ l in 
eq. iA5i . We assume efficient convection continues to just 
below the photosphere. 

At the photosphere, we adopt the usual Stefan-Boltzmann 
relation, F = aT 4 , and the hydrostatic condition, pn/A = 2/3, 
where A is the total acceleration defined in §|4] and 2/3 is the 
photospheric optical depth. Taking the photosphere to define 
the stellar surface, we compute the Lagrangian perturbations, 



and 



af _ 4 Ar 



Ap AA An 
p g K 



(30) 



(31) 



Using s and p as our independent thermodynamic variables, 
we write Ak/k = n^Ap / p + k s As jC p and AT /T = As/C p + 
V a[ jAp//7. In our numerical work (see § 15.4b . we self- 
consistently compute the perturbation AA to the effective sur- 
face gravity, in order to follow resonant oscillations, where the 
equilibrium-tide result fails. However, we are now addressing 
non-resonant forcing, for which we use the equilibrium-tide 
approximation at the surface, giving AA/g = -d^ r jdr (see 
§|4]l. We now have 



(32) 



(33) 



Ap _ ( K s As/C p + d£ r /dr 
P \ 1 + K ad 

and upon substitution, 

AF _ i f 1+K ad-VadK,A As 4V ad d& 
F \ 1+Kad / C p 1 + K ad dr 

Equation (g3) differs from iGoldreich & Wul d!999al) in that 
we retain the gravity perturbation in eq. ( TIT] ), whereas they 
consider a constant-gravity atmosphere (and no tidal pertur- 
bation). For g-modes in white dwarfs, the interesting region 
is near the surface and the motion is mainly horizontal, so that 
AA = Ag = is a good approximation. Since the equilibrium 
tide has large vertical motions, the AA term must be retained. 

The luminosity change across the convection zone is de- 
rived from the entropy equation (eq. I1A6II ). If we ignore 
horizontal flux perturbations (set I = in eq. IA6II ) and en- 
ergy generation, the equation for the luminosity perturbation 
AL/L = 2£ r /r+AF/F is 



d(AL/L) 
dM r 



= iuoTAs/L . 



(34) 



Integrating over the convection zone with constant As, we 
obtain 

AL ph 



AL h 



= iuAs / dM r T/L , (35) 

J cz 



where the subscript "ph" refers to the photosphere. We de- 
fine ? cz = Cpph J cz dM r T /L to be the mean thermal time of 
the convection zone, so that the right-hand side of eq. d35l > 
is iiot cz As/C P:Pb . 



Figure|5]shows that the thermal time at the base of the con- 
vection zone (of order f cz ) forM < 1 .3M is orders of magni- 
tude longer than the forcing periods of 1-10 days. Insofar as 
\AL\/L ~ |£,-|/r at any location in the star (i.e., if resonances 
are neglected), we see that (|£r|/r) _1 \As\/C p ~ (wf cz ) _1 < 1 
in stars with deep convective envelopes. In this limit, eq. d33l 
becomes 

AF 4V ad dtj r 

~E~~~TZ ~7T ■ ( 36 > 

r 1 + K a d or 

If we had set AA = 0, the amplitude of the photospheric flux 
perturbation would have been ~| As|/C p rather than the much 
larger value ~ | £ r \ /R . 

Photospheric flux perturbations in tidally forced solar- 
type stars with thick convective envelopes arise mainly from 
changes in the local effective gravity. This statement is remi- 
niscent of, but physically distinct from, von Zeipel's theorem 
(eqs. l2T1 and l22l ). We have recovered our equilibrium- tide 
scaling, AFg/F = —X^e./R, where eq. d36l l gives 

V ad 



Xt =4(1+2)- 



(37) 



1 + K ad 

For M = l.O-1.4M , we find A 2 ~ 1.9-1.1. These estimates 
neglect resonant excitation of g-modes, a point addressed in 

§m 

5.3. Analytic Result for Radiative Envelopes 

As the stellar mass increases beyond 1.4M Q , the outer 
convective region thins and sits close to the surface, where 
fed ~ fth- Figure [3] shows that the \ .6Mq model has two thin, 
inefficient surface convection zones, as well as a convective 
core. Radiative energy transport is important throughout the 
envelopes of these more massive stars. We now consider the 
idealized case of a completely radiative envelope, and obtain 
an analytic approximation for AL/L at the surface. 

Near the surface of a radiative star, we have H p /r <C 1, 
AuPp/M,- < 1, and u; 2 r/g < 1 for 2t:/lu = l-10days. Under 
these condition s, the quasi-adi abatic luminosity perturbation 
becomes (e.g. JUnno et al.lfT989h 



AL 



qad 



^Ap | g k\H p /y.d ( 



P H n 



where 



C = K ad-4V a d+- 



d\7 



ad 



(38) 



(39) 



V c/lnT ' 

and £,.. eq is the equilibrium-tide radial displacement (eq. lfT2l ). 
Nonzero values of Ap/p and (£, r - £,r,eq) / Hp indicate devia- 
tions from hydrostatic equilibrium. Care must be taken with 
these terms, because the denominators p and H p become very 
small close to the surface. 
With the help of the Appendix, we define the variables 

HpAp 

^ 

r p 



a=yi-yi+y3 = - 



(40) 



U H p Ap Sr-Zr, 

P=yi + — = — + 

gr r p r 



eq 



which satisfy the differential equations 
da d\np i gkf 
~dr~~~ 
d£ _ 
dr 



dr 

N 2 



uj- 



a+— . 
r 



(41) 



(42) 



(43) 
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FIG. 6. — Ratio of surface Lagrangian flux perturbation AF/F to 
equilibrium-tide displacement -U /gR for a range of forcing periods in the 
limit where surface g-modes are evanescent. The flux is evaluated at the loca- 
tion where uit^ = 1. Dashed, solid, and dotted curves correspond to M = 1.5, 
1.6, and 1.7Mq, respectively. 



When uj 2 <C gk 2 h H p , these equations produce the g-mode 
dispersion relation k 2 = k\N 2 /uj 2 (in the limit k 2 /k\ <C 1) 
for radial wavenumber k r . For these propagating waves, 
the surface amplitudes of a and j3 are determined at the 
core radiative-convective boundary , where g-modes are driven 
(e.g jGoldreich & Nicholson 1989). On the other hand, when 



^> gk\H p , the g-modes are evanescent (see lUnno et all 



1989) and we neglect the term gk\(3/Lo 2 in eq. 1421 . This 
limit yields the approximate solution a ~ —(£+4)(H p U /gR 2 ), 
or Ap/p ~ (4 + £)U / gR. In this case, Ap/p is not small 
compared to the fractional fluid displacement, and thus the 
equilibrium-tide approximation loses validity. 

From our stellar models, we find that the evanescent regime 
corresponds to forcing periods of <4-8 days for M = 1- 
1.6M Q , most of the range of interest. The high-frequency 
limit of eq. (l38l is 



L gR 



(44) 



This relation should be evaluated at the layer where Lut t h ~ 1, 
above which the luminosity effectively "freezes out." Fig- 
ure [6] shows the quasi-adiabic flux perturbation AF qa( j/F = 
AL qa( j/L-2^ r /R, evaluated where cofa = 1, for a range of forc- 
ing periods and M = 1.5-1.7M Q . Note that |AF/F| can be an 
order of magnitude larger than \U\/gR, because of the rather 
large values of \(\(£ + 4) for £ > 2. Much larger perturbations 
are possible when g-modes are resonantly excited in a radia- 
tive star, as we discuss in the next section. 

We must point out that the quasi-adiabatic approximation is 
technically inappropriate when ujt t ^ ~ 1 . Equation ( |44T > should 
be viewed as an estimate of the modulus of the luminosity per- 
turbation at the surface. If, for instance, |As|/C p > \U\/gR 
where ujt^ ~ 1, then AL/L at the surface will have a substan- 
tial imaginary part (see eq. Il34l ). This is what we find in the 
numerical calculations summarized in the next section. 



5.4. Numerical Examples 

Here we show solutions of the perturbed mass, momentum, 
and energy equations that describe linear, nonadiabatic oscil- 
lations of a star subject to a varying tidal f orce. The equations 
listed in the Appendix are the same as in lUnno et al.1 (fl989l) 
for radiative regions, but augmented to include the tidal ac- 
celeration. In convection zones, we apply the prescription 
As = constant based on our conclusions in § 15.11 Figure UJ 
summarizes how the interiors of 1M Q and 1.6M stars re- 
spond to resonant and non-resonant tidal forcing. The tidal 
potential has been scaled so that (, r /R = 1 corresponds to the 
equilibrium-tide surface displacement. 

For our 1 M Q model, the non-resonant response to a forc- 
ing period of ~3 days is shown in Fig. UJt. We see that £, r /R 
matches the equilibrium-tide result at the surface; the imag- 
inary piece is completely negligible. We also find that our 
approximation for AF/F at the surface (eq. l36l ) works very 
well. A factor of ^10 decay in \AL\/L occurred in order for 
| AF\/F ~ £ r /R at the surface. Variation of As/C p in the con- 
vection zone (\og[p I (GM 2 / R 4 )] > -2.5) is due to changes in 
C p . In the radiative interior, the oscillations are caused by 
most nearly resonant g-modes, whose amplitudes rise rapidly 
as the core is approached, due to conservation of wave lumi- 
nosity. We have checked that the quasi-adiabatic approxima- 
tion of AL/L is valid in the radiative region; the ratio of the 
real and imaginary parts is found t o be roug hly constant for 
the ingoing gravity-wave (see also Zahn 1975). 

In order to model the resonant response of a 1 M Q star, we 
tuned the forcing period to ~1 day (see Figs.UJ) and[8]). At the 
surface, both £ r and AL have dominant imaginary parts, due 
to the short radial wavelength of the g-mode compared to the 
equilibrium-tide fluid displacement. The entropy at the base 
of the convection zone is very strongly perturbed in compari- 
son to the non-resonant case, but AL is still damped by orders 
of magnitude as the surface is approached. 

Figure [8] shows the surface values of the complex modu- 
lus and phase of £ r /R and AF/F versus forcing period. The 
phase is tan" 1 (Imaginary /Real) S (— tt,tv). Solid lines connect 
points halfway between g-mode resonant periods. We find 
that the equilibrium-tide approximation given by eqs. (fl2l i and 
d36l l is excellent for non-resonant forcing. Dashed curves give 
the maximum and minimum values that occur on resonance. 
One example of a resonance is shown in the insets. Resonant 
forcing at periods of <2 days yields surface values of £, r /R 
and AF/F that differ substantially from the equilibrium-tide 
results. However, the ratio of resonance width to the spac- 
ing between adjacent resonances is ~10~ 4 , making resonant 
forcing very unlikely. It is noteworthy that at forcing periods 
of >2 days, the equilibrium-tide result holds extremely well 
even when precisely on a resonance. As explained by Zahn 
(1975), the resonant response can be considered as the sum of 
the equilibrium tide and the most nearly resonant wave. As 
the period increases, the g-mode radial wavelength decreases, 
resulting in a reduction of the overlap integral for the mode 
and the tidal force, which in turn gives a decreased amplitude 
of the wave component relative to the equilibrium tide. 

The non-resonant response of the 1.6M Q star is shown 
in Fig. UJ;. We see that the equilibrium- tide result pro- 
vides a good match to £, r /R- Our estimate for the modu- 
lus of the radiative luminosity perturbation in the evanes- 
cent limit (eq. fi4l ) agrees reasonably well with what is in 
Fig. UJ;. We also see that AL/L does roughly "freeze-out" 
when ojt t h ~ 1, just below the base of the convection zone at 
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FIG. 7. — Responses of tidally forced IMq and 1.6Mq main-sequence stars. Black, red, blue, and green curves denote, respectively, the logarithms of £ r /r, 
(6p + pSip)/ pgr, As/C p , and AL/L. Solid (dashed) curves show the real (imaginary) parts. The tidal potential has been scaled so that (,, /R = 1 corresponds 
to the equilibrium-tide value. The four panels show the following: (a) non-resonant response of a IMq star tidally forced at a period of 2tt/lu ~ 2.91 day; (b) 
resonant response of a IMq star with 2-n/u) ~ l.OOday; (c) non-resonant response of a 1.6Mq star with 2tt/ui ~ 3.00day; (d) resonant response of a 1.6M© 
star with 2tt /ui ~ 1 .02 day 

suits here, similar plots for masses between 1 M Q and 1 .6M Q 
show progressively more structure as the mass increases. The 
cause of this irregularity is not clear, but may have to do with 
the two thin surface convection zones changing the overlap of 
successive g-modes with the tidal force. 

6. SUMMARY 



log p J (GM 2 1 R 4 ) ~ -9 (see Pig . . Our expectations in § 
regarding the imaginary part of AL/L are borne out in Fig. [7]; 

A resonantly excited 1.6M Q star exhibits huge surface flux 
perturbations, radial displacements, and phase lags, as seen in 
Fig- EJi- In Fig. [9] surface values of \£, r \/R and \AF\/F are 
plotted as a function of forcing period, where we have taken 
care to resolve resonances. Resonant amplitudes vary non- 
monotonic ally with period, in contrast to the smooth behavior 
of the IMq star (Fig. |8). Although we do not show the re- 



We have investigated in detail the ellipsoidal oscillations of 
0.9-1.6M Q main-sequence stars induced by substellar com- 
panions. Classical models of ellipsoidal variability (e.g., 



10 



PFAHL, ARRAS, & PAXTON 



15 



10 




1.0S193 I 05193 I 05194 1. 05195 





forcing period [day] 



FIG. 8. — Surface radial displacement and Lagrangian flux perturbation 
versus forcing period for a 1 .QMq star. Solid lines connect points halfway be- 
tween resonant g-mode periods, while dashed curves give the maximum and 
minimum values found on resonance. The equilibrium-tide approximation is 
extremely good, except when the forcing period is <2 days and resonant. 



IWilsonlfl994l) are built on the assumption of hydrostatic bal- 
ance in a frame corotating with the binary orbit. This ap- 
proach is justified in the context of short-period (f or b < 
10 days) binaries containing two stars of comparable mass, 
where tidal dissipation circularizes the orbits and synchro- 
nizes the stellar spins with the orbital frequency. However, 
when the companion has a very low mass, we cannot assume 
that the binary is in complete tidal equilibrium; in fact, this 
state may be unattainable (see §[2]l. In this case, one must, in 
general, appeal to a dynamical description of the tidal inter- 
action. A substellar companion with P 01 -b > 1 day raises tides 
on the star that are a small fraction of the stellar radius (see 
e Q- QX permitting a linear analysis of the stellar oscillations. 

While the root of our study is a dynamical treatment of 
stellar tidal perturbations, the equilibrium-tide approximation 
does have an important realm of validity (see below). For 
this reason, we derived in §|3]a general expression (eq. fl6l ) 
for the measurable flux variation of a star that remains in 
hydrostatic equilibrium under the influence of a small exter- 
nal tidal force. This formula (1) assumes that the local per- 
turbation to the energy flux at the stellar surface is propor- 
tional to and in phase with the equilibrium-tide radial fluid 
displacement at each angular order £ (eq. Ifl2l ). (2) neglects 
stellar rotation, and (3) applies to inclined and eccentric or- 
bits. As expected, the fractional amplitude of the modulation 
is = (M p /M)(R/a) 3 for small eccentricities and/ = 90°, or 
~10" 5 (Mp/M y )(P orb /l day)" 2 for a star like the Sun (see §|2). 

A common practice is to use von Zeipel's theorem when 
computing the surface radiative flux from a tidally distorted 
star (see §|4j. The theorem assumes that the star is in hydro- 
static equilibrium and that the energy transport in the outer 
layers is purely by radiative diffusion. As already mentioned, 
the hydrostatic assumption is technically unjustified for sub- 
stellar perturbers. Moreover, the majority of Kepler targets 
will be main-sequence stars with masses of <1.4M Q , which 
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FIG. 9. — Surface radial displacement and Lagrangian flux perturbation 
versus forcing period for a 1.6Mq star. Curves connect evenly spaced points 
away from resonances, with finer spacing near resonance periods. 

have substantial surface convection zones. Evidently, von 
Zeipel's theorem is an inappropriate starting point for the con- 
ditions of inter est. 

Section |5TI discusses heat transport in perturbed stars with 
convective envelopes. Heuristic arguments are used to de- 
velop a simple treatment of the perturbed convection zone in 
main-sequence stars of mass < 1 .6M Q with forcing periods of 
1-10 days. We suggest that both the specific entropy s and its 
Lagrangian perturbation As are spatially constant i n convec- 
tive regions, a model partly inspired by the ideas of Brickhill 
(119831119901) . 

Using this prescription, we analytically compute in § l5.2l the 
perturbed flux at the photosphere of deeply convective stars 
(M < 1 AMq), where the thermal time scale at the base of the 
convection zone is much longer than the forcing period. We 
find that As/C p is negligible near the top of the convection 
zone, and that the photospheric flux perturbation is propor- 
tional to changes in the effective surface gravity. Thus, we re- 
cover the equilibrium-tide result, AF /F = -Xg^/R, at the sur- 
face, where depends on the adiabatic derivatives of opacity 
and temperature with respect to pressure (see eq. l37l ). Nu- 
merical solutions of the equations of linear, nonadiabatic stel- 
lar oscillations (see § !5.4l and Fig.|7^) corroborate our analytic 
estimates in the non-resonant regime. Resonant excitations of 
g-modes in the radiative stellar interior cause large departures 
from the equilibrium-tide approximation when the forcing pe- 
riod is <2days (Figs. [TJ^ and[8]l. However, the likelihood of 
being on a resonance is small, and at periods of >2 day the 
equilibrium-tide result holds for M ~ 1 Mq even with reso- 
nant forcing. 

Stars of mass >1.4M© have thin, relatively inefficient sur- 
face convection zones. Thus, g-modes can propagate very 
close to the surface and produce large flux perturbations and 
fluid displacements. Analytic arguments in § !5.3l indicate that 
the surface flux perturbations in these stars have non-resonant 
amplitudes of ^10e (eq. l44l and Fig. |6), in rough agree- 
ment with our numerical calculations (Fig. U). As seen in 
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TABLE 2 
Periodic Flux Modulations 



Variability Dominant Phase at 

Source Amplitude a,b Harmonic Maximum/Minimum c References 

Ellipsoidal 11 2 X lfT 5 nipmP^ 2 sin 2 / P orb /2 0.25(0.75)/0.00(0.50) 

Dopplei-e 3 x l(r 6 m p m-' 1 l 3 P^ 12 sin/ P olb 0.25/0.75 1 

Reflection f 6 X 10- 5 (a/0.1)m- 2 / 3 P~ 4/3 sin/ P orb 0.50/0.00 2,3 

Transit 10~ 2 m" 2 P olb ■■■ /0.00 4 



REFERENCES. — (1) Loeb & Gaudi 2003; (2) Seager, Whitney, & Sasselov 2000; (3) Sudarsky, Burrows, 
& Pinto 2000; (4) Seager & Mallen-Ornelas 2003 

a We assume that the orbit is circular in our estimates. b The dimensionless variables used are m p = 
M ( ,/(10~ 3 Mq), m = M/Mq and Pi = P mb /lday. We have assumed that the star and companion have re- 
spective radii of R/Rq = m and 0.1Pq. c The phase is in the range 0-1, where at phase the planet is closest 
to the observer. 1 ' Only the t = 2 component of eq. )16K with A2 = 2, is considered here. c We approximate 
the amplitude as 4v r /c, where v r is the reflex speed of the star along the line of sight, and the factor of 4 is 
approximately what one obtains for a V-band spectrum similar to the Sun. f Here a is the geometric albedo of 
the companion. The inclination dependence is an approximation for / ~ 90° and the Lambert phase function. 



Figs. [TJl and [9] a resonantly forced 1.6M Q star can exhibit 
flux perturbation amplitudes of >100e at forcing periods of 
~ 1 day. While the amplitudes are not as extreme at longer 
periods, their dependence on period is rather erratic (Fig. |9), 
an issue that deserves further study. It will be difficult to de- 
rive physical interpretations from the ellipsoidal variability of 
these more massive stars. 

7. DETECTION PROSPECTS 

The dominant sources of periodic variability of a star with a 
substellar companion are transit occultations (when cos/ < 
[R + R p ]/a), Doppler flux modulations, reflection of starlight 
from the companion, and ellipsoidal oscillations. For each 
of these signals, Table [2] lists the characteristic amplitude, 
period with the largest power in the Fourier spectrum, and or- 
bital phase(s) at which the light is a maximum or minimum. 
The transit contribution is included for completeness, but its 
duration is sufficiently short — a fraction ~{R + R p )/(tto) of 
P 01 b — th at it should often be pos sible to excise it from the 
data (see lSirko & Paczvris ki 2003). Of the remaining signals, 
the Doppler variability is the simplest, being purely sinusoidal 
with period P 01 -b when the orbit is circular. The dominant I = 2 
piece of the equilibrium-tide approximation to the ellipsoidal 
variability (see eqs. |[T6l and fTTl ) is also sinusoidal when 
e = 0, but with period P 01 -b/2. Reflection is more problem- 
atic, as its time dependence is generally not sinusoidal and 
not known a priori. 

If the companio n scatters light as a Lambert sphere (e.g., 
Seage r~et~aTI l2000b . the Fourier spectrum of the reflection 
variability has finite amplitude at all harmonics of the orbital 
frequency il, but the amplitude at 2il is roughly 1/5 of the 
amplitude at fi. Therefore, the reflection and ellipsoidal vari- 
ability amplitudes may be similar at a frequency of 20 when 
a = 0.1, M p ~ My, and P or b ~ 1 day. Also, the orbital phase 
at which the reflected light is a maximum is distinct from 
both the Doppler and ellipsoidal cases, further distinguishing 
the signals. However, Lambert scattering is probably never 
appropriate in real planetary atmospheres. Infrared reemis- 
sion of absorbed optical light, multiple photon scattering, and 
anisotropic scattering typically conspire to narrow the peak in 
the reflection lightcurve and lower the albedo, decreasing the 
prominence of the reflection signal. These issues are sensi- 
tive to the atmospheric chemistry and the uncertain details in 
models of irradiated giant planets. For reasonable choices re- 



garding the atmospheric composition, calculated o ptical albe- 
dos of Jovian planets range from <0.01 to ~0.5 dSeager et al.l 
2000; Sud arsky et al.l2000h . Recent photometric observations 
of HD 209458, the star hosting the first-detected transiting gi- 
ant planet ( P OI h, ~ 3.5 days), constrain the planetary albedo to 
be <0.25 dRowe et al.ll2006l) . 

Detailed lightcurve simulations will be required to say how 
well the different periodic signals can be extracted from the 
data. This is beyond the scope of the current study. We now do 
the simpler exercise of isolating the ellipsoidal modulations 
and assessing when this effect alone should be detectable. For 
a star of apparent visual magnitude V and an integration time 
of T = 6hr, Kepler's photon shot noise is 7 

Instrumental noi se should contribute at a similar level (e.g., 
iKoch et al.ll2006l) . If the data is folded at the orbital period 
and binned in time intervals T <C P or b, the shot noise is sup- 

-1/2 

pressed by a factor of ~« orb , where n 01 -b is the number of 
folded cycles. After folding 1 year of continuous photomet- 
ric data using T = 6hr, a star with V < 12 orbited by a giant 
planet with f or b < 3 days may have a fractional shot noise per 
time bin of <10~ 6 . This is less than the ellipsoidal amplitude, 
((5J?"/J£")eii, when / is not too small. 

The actual situation is not so simple when the data spans of 
weeks or months, because the intrinsic stochastic variability 
of the star will not have a white-noise power spectrum. Over 
times of <1 day, the Sun shows variability of (<5J?/J^")i n t ~ 
10~ 5 , but the amplitude rises steeply between ~1 and 10 days 
to ^10~ 3 . Intrinsic variability tends to be large near the rota- 
tion period of the star, due mainly to starspots. Low-frequency 
variability may not too damaging for the study of ellipsoidal 
oscillations induced by planets with P or b < 3 days, but more 
study is needed. 

Kepler's target list will contain ~10 5 main-sequence FGK 
stars with V = 8-14. The statistics of known exoplanets in- 
dicate that 1-2% of all such st ars host a giant pl anet (M p > 
Mj) with P orb < lOdays (e.g.. iMarcy et all [2005). Of these 
"hot Jupiters," ~30% have P or b = 1-3 days. It seems that a 

7 An integration time of T = 6 hr is chosen for convenience; Kepler's nom- 
inal exposure time is 30 min. Here we use the V-band flux as a reference, but, 
in fact, the Kepler bandpass is 430-890 nm, which spans B, V, and R colors. 
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TABLE 3 

Number of Kepler Stars with Detectable 
Ellipsoidal Oscillations 
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12 
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33 


19 


7 
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1 



maximum of ^10 3 Kepler stars will have detectable ellip- 
soidal modulations. If we neglect intrinsic stellar variabil- 
ity and consider only shot noise, then many systems with 
Ptxb ^ 3 days and V < 14 will have signal-to-noise S/N > 1 
after ~100 cycles are monitored; this may amount to >100 
stars. Obviously, the number drops when we place higher de- 
mands on S/N and include the intrinsic variability. The results 
depend critically on the distributions of M p and P rb- 

In order to better estimate the number of stars with poten- 
tially detectable ellipsoidal oscillation, we perform a simple 
population synthesis calculation. Denote the set of star-planet 
system parameters by P = {M,M p ,P m ±,,I}, and let /(P)c/P be 
the probability of having a system in the 4-dimensional vol- 
ume dP. We assume that the planetary orbits are circular and 
obtain (5jF/J^") e ii from the equilibrium-tide estimate in Ta- 
ble |2] Given the mass of the star, we compute its absolute 
V magni tude on the main-sequenc e using the approximation 
(see also lHenry & McCarthyll 19931) 



= 4.8- 10.3 log(M/M ) , 



(46) 



which is in accord with the usual mass-luminosity relation 
log(L/L©) ~ 41og(M/M Q ) for M ~ 1 M Q . With a maximum 
apparent magnitude of V max = 14 for the Kepler targets, the 
maximum distance of the star is 



Anax = 10 



1+0.2(14-^1,) 



pc 



(47) 



With a certain signal-to-noise threshold (S/N) m \ n , there is a 
maximum distance D d < D, mx to which the ellipsoidal vari- 
ability is detectable. For a spatially uniform population, the 
detectable fraction of systems is (Dd/D max ) 3 . Thus, the net 
detectable fraction among all systems is 



S= / dPf(P) 



Dd 



(48) 



an integral over all relevant P space. 

When the only noise is intrinsic to the star, N = (<5J^/J^)i nt 
and S/N is independent of distance, so that Dj ; nt = £> max when 
S/N > (S/AOmin, an d A/.int = otherwise. In the case of pure 
shot noise, there is a maximum magnitude for which the 
ellipsoidal oscillations are detectable: 



V d = 51og 



ell 



X (S/N) n 



(49) 



where \ ~ 10~ SA (T 6 n m T 1/2 is the value of (<^/J0shot for 
V = 0, T = 67f,hr, and « rb = lOOnioo. The corresponding dis- 
tance is given by log[A/. s h o t/10pc] = 0.2(Vrf — „#y) if Vd < 
VW, and is D d . sho( = D max when V d > V max . We take the max- 
imum detectable distance to be D d = mm{Dd ,i n t, A/ .shot}- 

At this point the simplest approach is to assume that the 
parameters {M,M p ,P OT b,I} are statistically independent and 



carry out a Monte Carlo i ntegration to obtain S. To this 
end, we draw M from the Kroupa et al] ( 1 19931) initial mass 
function in the range of 0.5-1.5Mq. The planetary mass is 
chosen from the di stribution f(M p ) oc M~ x for M p = l-10My. 
Marcy ~e7aT1 (120051) find that x ~ 1 when considering all de- 
tected planets; the shape of f(M p ) is not well constrained at 
P orb < 10 days. We let x = 1 and 2. We adopt f(P mh ) oc P~~ h 
over 1-10 days. Multiplying the resulting value of S by 1000 
provides a crude estimate of the actual number of Kepler tar- 
gets with detectable ellipsoidal variability. No single value of 
y is consistent with the data, and so we consider the reason- 
able range y = — 1, 0, and +1. Inclinations are chosen under 
the assumption that the orbits are randomly oriented, such that 
/(cos/) =1/2 for / £ (0,7r). Our calculations use fixed values 

of (<JJ7JDtat = 10" 5 and T b = n 100 = 1. 

Results of our Monte Carlo integrations are shown in Ta- 
ble|3]as actual numbers of Kepler targets. The largest number 
of detectable systems is obtained when x = y = 1, parameters 
that yield the largest proportions short periods and massive 
planets. We expect that ^10-100 Kepler stars may exhibit 
ellipsoidal oscillations with S/N > 5. A handful of systems 
might have S/N > 10. Higher harmonics from the I = 3 com- 
ponents of eq. (fTot or modest eccentricities might be accessi- 
ble for at most a few stars. 

Our integrations also check for cases where the planet is 
transiting. As (S/N) m i„ increases from 1 to 5, the fraction 
of systems in Table [3] with |cos/| < (R + R p )/a runs from 
~30% to ~50%, with a weak dependence on x and y. Such 
significant fractions stand to reason, since systems with the 
shortest periods have the highest ellipsoidal amplitudes and 
transit probabilities. Transit measurements directly give f or b, 
sin/ > 0.95 (for P orb > 1 day), and (R p /R) 2 . The planet mass 
M p can be determined with the addition of spectroscopic ra- 
dial velocity measurements, which should be possible for 
most of the Kepler targets with detectable ellipsoidal oscil- 
lations. The ellipsoidal amplitude then depends on the un- 
measured stellar mass and radius via e oc R 1 /M 2 (eq. JT]), as 
well as the stellar photospheric conditions (eq. 1361 ). If M 
and R are obtained from stellar models, ellipsoidal variabil- 
ity may provide an interesting consistency check on all the 
system parameters, as well as test the theory of forced stellar 
oscillations. 

As a last point, we emphasize that stars of mass >1.4M Q 
may have typical ellipsoidal amplitudes of ~10e. However, 
such stars will also be younger than most Kepler targets and 
probably have intrinsic variability ^10~ 5 . We carried out 
Monte Carlo integrations with M = 1 .4-1 .6M Q , (5JP/JF) e n = 
lOesin 2 /, and x = y = \. As we vary /^)wt from 10~ 5 
to 10~ 4 , $ decreases from large values of ~0.4 to a small 
fraction of ~0.03 for (S/A0min = 10. Unfortunately, we do 
not know how many such stars will be included in the Kepler 
target list. Also, there has not yet been a discovery of a giant 
planet with f rb < 10 days around a star of mass > 1 AMq, but 
exoplanet surveys tend to exclude these more massive stars. 



We thank Tim Brown for general discussions and address- 
ing Kepler questions, J0rgen Christensen-Dalsgaard for guid- 
ance on stellar luminosity perturbations, and Mike Muno for 
advice on signal processing. This work was supported by NSF 
grant PHY05-5 1164. 



ELLIPSOIDAL OSCILLATIONS 



13 



APPENDIX 

OSCILLATION EQUATIONS 

Her e we list the nonadiabatic, linearized fluid equations that we solve numerically. The reader is referred to lUnno et all 
(1989) for a complete discussion. Scalar and vector quantities are expanded in spherical harmonics Y(, m and poloidal vector 
harmonics, respectively. The momentum, mass, and energy equations are written in terms of the dimensionless variables y\ = £ r / r, 
y 2 = (5p/p + 5tp)/gr, j3 = Sip/gr, y^ = g~ x d8ip/dr, y$ = As/C p , and ys = AL/L. Here L is the total (radiative plus convective) 
luminosity. The radial flux perturbation is AF/F = AL/ L-2^ r /r. In radiative zones, the nonadiabatic equations are 



dy\ 
dlnr 

dyi 
dlnr 

dyi 
dlnr 

dyA 
dlnr 

dys 
dlnr 



--y\ 



8 --3 

cl 



+yi 



•5 

cf / cf UJ A 



-N 2 



g/r 

: 3 ; 3(l-?y)+3'4 , 
N 2 gr 

: yiv—r+y2r]—+y3 

g r cf 



+y 2 1-77+- 



AT 



N 



/rj 8/ 



-p s y5- 



l du 

g dr 



=y\ 



H, 



V ad T) 



+ 4(V-V ad ) + C 2 



+yi 



(Vad-V) 



8 k l r 



■c 2 



r r r r r 

+B— C2+J4— Vad+ys— V(4-K J )-y 6 — V+ — 
tip tip tip tip tip 



fdU/dr k\ 



U 



dye 
dlnr 



Vad 

V 



-yiW+D^+yiW+D^+ys 



4irr 3 pC p T £(£+l)H p 



(Al) 
(A2) 
(A3) 
(A4) 

(A5) 
(A6) 



where c s is the sound speed, r] = dlnM r /dlnr, c 2 = (r/// ;; )V(K a d _ 4V a d)+ V a d(^ln V a i/dlnr+r/H p ), and we have ignored energy 
generation terms. Note that the tidal acceleration has been added to the momentum equations. In convection zones, we 

ignore turbulent viscosity effects and replace the radiative diffusion equation (eq. IA5I ) with the prescription As = constant (see 
§ 15. It . or more precisely 



(A7) 



Equation ( IA6b still involves the total (convective plus radiative) luminosity. We ignore energy generation and horizontal flux 
perturbation terms, i.e. we ignore all terms with spherical harmonic index £ in eq. (1A61 > in convection zones. 

At the center of the star, we require the solutions to be finite, and also set As = 0. At the surface, we set dp = pg^ r and we 
require Sip to decrease outward. This boundary condition is only approximate, as g-modes may propagate above the convection 
zone for wave periods of >4days in our 1M Q model. The final surface boundary condition is given by eq. OTb . Care must be 
used in the radiative zone just below the photosphere, since the entropy perturbation is far from the quasi-adiabatic value. If we 
solve the radiative diffusion equation in this region, we find that the entropy increases by ~10 orders of magnitude in just a few 
grid points. However, we regard this behavior as unphysical, because the region at the top of the convection zone is optically thin. 
To eliminate this unphysical behavior, we set As to a constant at such low optical depths. 
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